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Abstract 

The Lagrangian theory of gravitational instability of homogeneous-isotropic Friedman-Lemaitre 
cosmogonies investigated and solved in the series of papers by Buchert (1989), (1992), Buchert & 
Ehlers (1993), Buchert (1993a, b), Ehlers & Buchert (1993), is illustrated. The third-order solution 
of this theory for generic initial conditions is presented and realized in a special case by employing 
methods of high-spatial resolution of the density field. 



1 A class of generic third-order solutions 

We require that, initially, the peculiar-acceleration w{X,t) be proportional to the peculiar- velocity 
u{X,t): 

u{X,to) = w{X,to)to , (1) 

where we have defined the fields as usual (compare Peebles 1980, Buchert 1992). Henceforth, we 
denote the peculiar- velocity potential at the initial time tg by S, u(X,to) =■ VqS, and the peculiar- 
gravitational potential at tg by cj), w(X,io) =: — Vo<?!>, where Vq denotes the nabla operator with 
respect to Lagrangian coordinates X . This restriction has proved to be appropriate for the purpose 
of modeling large-scale structure, since, for irrotational flows, the peculiar- velocity fleld tends to be 
parallel to the peculiar-acceleration after some time, related to the existence of growing and decaying 
solutions in the linear regime. Still, any generic function (S or cj)) can be given initially. 

With a superposition ansatz for Lagrangian perturbations of a flat homogeneous and isotropic 
background the following family of trajectories x = f(X, a) as irrotational solution of the Euler-Poisson 
system up to the third order in the perturbations from homogeneity has been obtained (Buchert 1993b). 
The general set of initial conditions ((^(X), (5(X)) is restricted according to S = —cj) tg (see equation 
(1)). (The parameter e is thought to be absorbed into the amplitudes of the perturbation potentials; 
a{t) = {tltgfl^y. 



with: 



qzz 



3 
2 
3 
2 
3 
2 
3 
2 
3 
2 

where the initial displacement vectors have to be constructed by iteratively solving the 7 Poisson 
equations: 

Ao5(^) = 2/7(5^) 
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(i, J, A; = 1, 2, 3 cyclic ordering). 

The scalar potential iS'^^'^ and the vector potential S'--^'^^ generate interaction among the first- and 
second-order perturbations. Note that the general interaction term is not purely longitudinal. Inspite 
of the irrotationality of the fiow in Eulerian space, vorticity is generated in Lagrangian space starting 
at the third order. 

An important remark relevant to any realization of the solution (2) concerns the possibility of set- 
ting (S*^^) = Sto without loss of generality, if periodic boundary conditions are imposed on S (compare 
Buchert 1992, Ehlers & Buchert 1993). In the latter paper it is proved that the solution (2) is unique 
for such boundary conditions. With this setting, the first-order solution reduces to the well-known 
"Zel'dovich approximation" (1970) as the only model which depends locally on the initial conditions. 

The presented model is currently being compared with numerical simulations. We analyze cross- 
correlations for various power spectra with PM simulations (Buchert et al. 1993 a,b), and special 
cluster models at high resolution using a tree-code (Buchert et al. 1993 c). 

Parallel attempts to investigate Lagrangian perturbation solutions are made by Moutarde et al. 
(1991), Bouchet et al. (1992) and Lachieze-Rey (1993). 



2 Realization at high-spatial resolution 

Let us use for the illustration a coherent random fiuctuation field as initial condition. We specify this 
to standard HDM initial conditions in a box of size 50Mpc h~^ normalized to cig = 1 using top-hat 
smoothing. 

The solution (2) is realized for the first- and third-order approximations at this amplitude (see 
Figs.l). Since the collapse process is significantly accelerated by the higher-order corrections, the 
first-order approximation represents an earlier stage in the evolution of the clustering process. 




(a) (b) 

Figure 1: A thin slice through a 3D HDM realization at high-spatial resolution (1024^ trajectories 
collected into a 128^ pixel grid) is shown for the first-(a) and the third-order approximation (b) for 
the evolution stage z = for the normalization o-g = 1. 

The models presented in Fig.l have been realized at high-spatial resolution employing the methods 
described in Buchert & Bartelmann (1991). In particular, we calculate the initial condition at 128^ 
grid points and interpolate 12 times trilinearly between the grid points to obtain an effective resolution 
of 1024^ particles. This interpolation yields highly accurate results, since the initial fluctuation field 
is sufficiently smooth and the mapping is analytical providing formally infinite resolution. 
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